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ABSTRACT  Temperature  is  important  for  mosquito  development  and  physiological  response,  Sev¬ 
eral  genes  of  heat  shock  protein  (HSP)  families  are  known  to  be  expressed  in  mosquitoes  and  may 
be  crucial  in  responding  to  stress  induced  by  elevated  temperature.  Suppression  subtractive  hybrid¬ 
ization  (SSH)  was  used  to  identify  target  transcripts  to  heat  shock  treatment  in  female  Aedes  aegypti. 
Subtraction  was  performed  in  both  directions  enriching  for  cDNAs  differentially  expressed  between 
a  non-  heat  shock  control  and  heat  shock  treatment.  Heat  shock  treatment  of  female  Ae.  aegypti  was 
carried  out  for  1  h  at  42°C.  Clones  from  differentially  expressed  genes  were  evaluated  by  sequencing. 
Target  transcripts  up-regulated  by  heat  shock  included  five  different  HSP  gene  families  and  27  other 
genes,  such  as  cytochrome  c  oxidase,  serine-type  endopeptidase,  and  glutamyl  aminopeptidase. 
Additionally,  some  novel  genes,  cytoskeleton  and  ribosomal  genes,  were  found  to  be  differentially 
expressed,  and  three  novel  up-regulated  sequences  belonging  to  a  low-abundance  class  of  transcripts 
were  obtained.  Up-regulated /down-regulated  transcripts  from  heat  shock  treatment  were  further 
confirmed  and  quantified  by  quantitative  real-time  polymerase  chain  reaction  (PGR),  High  temper¬ 
atures  can  alter  the  gene  expression  of  a  vector  mosquito  population,  and  further  characterization  of 
these  differentially  expressed  genes  will  provide  information  useful  in  understanding  the  genetic 
response  to  heat  shock  treatment,  which  can  be  used  to  develop  novel  approaches  to  genetic  control, 
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Higher  temperatures  can  drastically  alter  the  genetic 
structure  and  gene  expressions  of  a  vector  mosquito 
population.  Several  families  of  heat  shock  proteins 
(HSPs)  are  known  to  be  expressed  in  insects  and  may 
have  a  cumulative  role  in  determining  stress  in  re¬ 
sponse  to  elevated  temperature  (Mahroof  et  al.  2005; 
Yadav  et  al.  2005;  Rinehart  et  al,  2006a,  b;  Robich  et  al, 
2007),  which  itself  is  governed  by  several  genes.  In 
addition  to  heat  shock,  the  expression  of  HSPs  has  also 
been  reported  to  be  induced  under  various  stress  con¬ 
ditions  (such  as  pathogen  infection,  heavy  metal  ions, 
hypoxia,  and  osmotic  stress)  in  many  other  animals 
and  insects  including  mosquitoes  (Mosser  et  al,  1988; 
Yamuna  et  al.  2000;  Boone  and  Vijayan  2002a,  b;  Spees 
et  al.  2002;  Gheng  et  al.  2003;  Sim  et  al.  2005;  Ghuang 
et  al.  2007;  Sim  et  al.  2007),  HSPs  are  also  develop- 
mentally  regulated.  For  example,  the  synthesis  of  HSP 
in  Anopheles  stephensi  was  correlated  with  the  various 
morphological  and  physiological  events  occurring 
during  development  (Gakhar  and  Shandilya  1999). 
Therefore,  HSPs  are  important  for  us  to  better  under¬ 
stand  the  mechanism  of  heat  shock  and  other  types  of 
stress  in  mosquitoes. 


^  Corresponding  author,  e-mail;  liming.zhao@ars.usda.gov. 


In  this  study,  SSH  was  used  to  identify  genes  that  are 
differentially  expressed  during  heat  shock  treatment 
of  female  Aedes  aegypti  L.  Expression  patterns  were 
confirmed  by  quantitative  real-time-polymerase  chain 
reaction  (RT-PGR),  and  the  regulated  genes  identi¬ 
fied  are  discussed  in  the  context  of  their  possible 
function  at  high  temperature  (heat  shock  treatment) . 
As  part  of  our  effort  to  develop  new  insecticides  for 
mosquito  control,  the  finding  of  HSP  genes  and  genes 
from  other  families  specifically  expressed  in  heat 
shock-treated  female  Ae.  aegypti  may  provide  infor¬ 
mation  needed  to  identify  proteins  critical  to  mosquito 
survival.  Using  RNAi  technology  to  knock  down  these 
critical  proteins  may  provide  additional  targets  that 
can  be  developed  as  molecular  pesticides  (Pridgeon  et 
al.  2008). 

Materials  and  Methods 

RNA  Extraction.  Aedes  aegypti  ( Orlando,  EL,  strain, 
maintained  since  1952)  were  reared  in  the  insectary  of 
the  Mosquito  and  Fly  Research  Unit  at  the  Genter  for 
Medical,  Agricultural,  and  Veterinary  Entomology, 
USDA-ARS,  Gainesville,  FL,  To  test  how  tempera¬ 
tures  affect  gene  expression  in  mosquitoes,  extreme 
heat  shock  treatment  (42°G)  of  female  Ae.  aegypti  was 
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conducted  to  optimize  the  response  and  facilitate  de¬ 
tection  of  the  highly  expressed  genes.  Specifically, 
10-d-old  females  were  exposed  for  1  h  at  42°C  and  56  ± 
1.5%  RH  in  an  environmental  chamber  (L-C  Incuba¬ 
tor;  Lab-Line  Instruments,  Melrose  Park,  IL)  for  this 
study.  Untreated  females  were  held  at  a  constant  room 
temperature  (23°C)  and  56  ±  1.5%  RH  in  an  environ¬ 
mental  chamber.  The  total  RNAs  were  extracted  using 
TRIzol  reagent  according  to  the  manufacturer’s  in¬ 
structions  (Invitrogen,  Carlsbad,  CA) .  Poly  (A)  ^  RNA 
was  isolated  applying  Oligotex-dT  suspension  (QIA- 
GEN,  Valencia,  CA) .  RNA  samples  were  quantified  by 
SmartSpec  Plus  Spectrophotometry  (HIO-RAD,  Her¬ 
cules,  CA). 

Subtractive  Hybridization.  Differentially  expressed 
mRNA  populations,  both  from  heat  shock  treatment 
and  from  untreated  controls,  were  converted  into 
cDNA  using  PCR-Select  cDNA  Subtraction  Kit  (Clon- 
tech;  Roche  Molecular  Systems,  Alameda,  CA).  Sub¬ 
tractive  cloning  is  a  powerful  technique  that  allows 
isolation  and  cloning  of  mRNAs  differentially  ex¬ 
pressed  in  two  different  populations.  In  the  general¬ 
ized  subtraction  scheme,  the  mosquitoes  to  be  com¬ 
pared  are  the  [  +  ]  or  tracer  (HS  treatment)  and  the 
[  — ]  or  driver  (non-HS,  control),  where  mRNAs  ex¬ 
pressed  in  the  tracer  and  not  the  driver  are  isolated. 
Hybrids  that  form  include  sequences  common  to  both 
mosquito  populations.  The  unhybridized  fraction  is 
enriched  for  sequences  that  are  preferentially  ex¬ 
pressed  in  the  tracer  mosquito  population.  Differen¬ 
tially  expressed  cDNA  sequences  are  used  to  con¬ 
struct  a  subtracted  cDNA  library  (Patel  and  Sive 
2001). 

PCR-select  cDNA  Subtraction  Library.  Forward 
and  reverse  subtracted  libraries  were  cloned  using  the 
TOPO  TA  Cloning  Kit  for  Sequencing  (Invitrogen). 
Transformed  plasmids  were  inserted  into  One  Shot 
TOPIO  Competent  Cells  (Invitrogen)  and  grown 
overnight  on  Luria-Rertani  (LE)  plates  containing 
ampicillin  and  X-Gal  (5-bromo-4-chloro-3-indolyl- 
beta-D-galactopyranoside).  For  each  library,  >100 
white  colonies  were  isolated  and  grown  overnight  in 
LB-ampicillin  broth  at  37°C  and  235RPM  in  the  Innova 
4000  Incubator  Shaker  (New  Brunswick  Scientific, 
Edison,  NJ). 

Gene  Sequencing  of  PCR-select  cDNA  Subtraction 
Library.  Clones  of  the  subtracted  library  were  purified 
with  QIAprep  Miniprep  (QIAGEN).  The  plasmid 
DNAs  were  digested  by  using  EcoRI  enzyme  for  1.5  h 
and  were  run  on  a  1%  agarose  gel  to  confirm  the  DNA 
insert.  Selected  clones  were  sent  to  the  DNA  Sequenc¬ 
ing  Core  at  the  Interdisciplinary  Center  for  Biotech¬ 
nology  Research  (ICBR),  University  of  Florida,  to  be 
sequenced,  and  the  sequences  were  analyzed  using 
the  National  Center  for  Biotechnology  Information 
(NCBI)  BLASTN  program  to  identify  sequence  ho¬ 
mologies.  The  homologous  DNA  fragments  in  the 
gene  or  cDNA  or  mRNA  are  also  recorded  in  Table  1. 

cDNA  Synthesis  for  Real-time  PCR.  A  5-p,g  aliquot 
of  purified  RNA  was  reverse  transcribed  in  20-fd  re¬ 
action  volume  using  Superscript  II  first-strand  cDNA 
Synthesis  system  for  RT-PCR  according  to  the  man¬ 


ufacturer’s  instruction  (Invitrogen) ,  The  reaction  was 
terminated  by  heat  inactivation  at  95°C  for  5  min.  The 
cDNA  samples  for  heat  shock  treatment  and  control 
were  diluted  by  adding  80  /d  ddHjO  (500  ng/  fA)  and 
stored  at  — 20°C. 

To  design  gene-specific  primers,  detailed  analyses 
of  the  nucleotide  sequence  of  all  up- /down-regulated 
genes  were  performed  using  PRIMER3-Design  Primer 
Pairs  and  Probes  program  from  Biology  Workbench 
(http://workbench.sdsc.edu) .  The  primers  for  ACTIN 
gene  were  also  designed  for  an  internal  control  and 
comparison.  The  primers  are  listed  in  Supplemental 
Table  1. 

Real-time  PCR  Amplification.  The  real-time  PCR 
assay  for  AeaGene  was  performed  using  Platinum 
SYBR  Green  qPCR  SuperMix-UDG  with  ROX  (In¬ 
vitrogen)  in  a  volume  of  15  /xl  on  an  Applied  Biosys¬ 
tems  7300  Fast  Real-Time  PCR  System  (Foster  City, 
CA) .  The  PCR  mixture  consisted  of  1  p-1  diluted  cDNA 
(.500  ng/ fjl),  0.5  p-M  primers,  and  IX  master  mix.  In 
every  RT-PCR  run,  ACTIN  was  used  as  an  internal 
control  to  normalize  for  variation  in  the  amount  of 
cDNA  template.  The  PCR  primers  used  were  Aea- 
GENE-F  and  Aea-GENE-R  (Supplemental  Table  2). 
The  PCR  primers  for  ACTIN  and  the  PCR  thermal 
cycling  parameters  were  the  same  as  described  in  the 
previous  publication  (Zhao  et  al,  2008).  This  experi¬ 
ment  was  replicated  three  times.  Relative  expression 
levels  were  calculated  as  follows.  First,  AeaGENE  tran¬ 
script  levels  relative  to  a  standard  (ACTIN)  using  the 
formula  AC-r  =  C^-  (AeaGENE)  -  C^-  (ACTIN) .  Second, 
an  average  AC^  value  for  each  sample  was  calculated. 
Third,  relative  expression  levels  were  calculated  using 
the  equation  100  X  '^‘^1 .  Last,  fold  increases 

or  decreases  were  calculated  using  HS  treatment  RNA 
relative  expression  level /non- heat  shock  treatment 
control  RNA  relative  expression  level. 

Results 

Identification  of  Genes  Specifically  Expressed  Dur¬ 
ing  Heat  Shock  Treatment  in  Ae.  aegypti.  To  show 
how  global  genes  are  specifically  expressed  during 
heat  shock  treatment  in  Ae.  aegypti,  we  used  subtrac¬ 
tive  hybridization  to  successfully  construct  a  subtrac¬ 
tion  cDNA  library  of  control  female  mosquitoes  and 
female  after  heat  shock  treatment  for  1  h  at  42°C. 
Clones  of  differentially  expressed  genes  were  evalu¬ 
ated  by  sequencing.  Using  NCBI  blast  to  analyze  the 
sequences,  we  found  10  clones  encoding  Ae.  aegypti 
heat  shock  protein  genes,  including  at  least  three 
groups:  HSP  26  kDa,  other  HSPs,  and  heat  shock  cog¬ 
nate  (HSC)  70  mRNA.  A  full-length  cDNA  clone 
encoding  an  HSP  26  kDa  family  was  isolated  and 
submitted  to  National  Center  for  Biotechnology 
Information  (NCBI  accession  no.  gb|EU048558.l|). 
There  were  12  clones  encoding  10  non-HSP  genes. 
They  were  Ae.  aegypti  serine- type  endopeptidase,  glu¬ 
tamyl  aminopeptidase,  cytochrome  c  oxidase  (subunit 
VIA),  casein  kinase,  adp/ATP  carrier  protein,  elon¬ 
gation  factor  2,  ATP  synthase  delta  chain,  thioredoxin 
reductase,  methionine  adenosyltransferase,  and  trans- 
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Table  1.  Heat  shock  treatment  response  genes  from  Ae.  aegypti  isolated  by  suppressive  subtractive  hybridization  with  homologous 
gene  blast  from  NCBI 


Clone 


Accession  no.,  organism,  putative  identity,  mRNA/cDNA  length 


Gene  regions 


HSPs 
51,  61 
62 
50 

12,  21 

41 
22 

44,  60 

Non-HSPs 

9,  26 

10,  17 
14 

16 

20 

36 

37 

42 
45 

49 

Unknown 

18 

25 

35,  59,  64,  68 

38 

39 

47,  52,  55,  70,  71 


58 

3,  32 

4,  29,  66 

Cytosolic  large 

ribosomal  subunit 
1,  69 

27 

43 

40S  ribosomal  protein 
2 
30 

67 

60S  ribosomal  protein 
5 
8 
11 


24 
33 
40 
Actin 
31,  34 


gb  EU048558.1,  Ae.  aegypti,  HSP  26-kDa  mRNA,  complete  cds,  length  =  738 
ref  XM_001655592.1,  Ae.  aegypti,  HSP  partial  mRNA,  length  =  2,605 
ref  XM_001649702.1,  Ae.  aegypti,  HSP  partial  mRNA,  length  =  2,602 
ref  XM_001649702.1,  Ae.  aegypti,  HSP  partial  mRNA,  length  =  2,602 

gb  AY432606.1,  Ae.  aegypti,  putative:  HSP  mRNA  sequence,  length  =  1,275 
gb  DQ440299.1,  Ae.  aegypti,  HSC  70  mRNA,  complete  cds,  length  =  1,956 
gb  DQ440299.1,  Ae.  aegypti,  HSC  70  mRNA,  complete  cds,  length  =  1,956 


ref  XM_001659911.1,  Ae.  aegypti,  serine-type  enodpeptidase  partial  mRNA.  length  = 
1,117 

ref  XM  OO 1658502.1,  Ae.  aegypti,  casein  kinase  partial  mRNA,  length  =  3,840 
ref  XM_001658180.1,  Ae.  aegypti,  glutamyl  aminopeptidase  partial  mRNA,  length  = 

3,060 

ref  XM_001656578,1,  Ae.  aegypti,  cytochrome  c  oxidase,  subunit  VIA  partial  mRNA, 
length  =  522 

ref  XM_001654983.1,  Ae.  aegypti  adp,  ATP  carrier  protein  partial  mRNA,  length  =  835 
gb  AF331798.1  AF331798,  Ae.  aegypti  elongation  factor  2  (Ef-2)  mRNA,  complete  cds, 
length  =  2,681 

ref  XM_001654955,1,  Ae.  aegypti,  ATP  synthase  delta  chain  partial  mRNA,  length  =  902 
ref  XM  OO 16626 16,1,  Ae.  aegypti,  thioredoxin  reductase  pai'tial  mRNA,  length  =  2,137 
gb  AY432324.1,  Ae.  aegypti,  putative:  methionine  adenosyltransferase  mRNA  sequence, 
length  =  2,655 

gb  AY431275.1,  Ae.  aegypti,  translational  inhibitor  mRNA  sequence,  length  =  663 

gb  AY432385.1,  Ae.  aegypti,  conserved  unknown  mRNA  sequence,  length  =  1,198 
ref  XM  OO 1660484,1,  Ae.  aegypti,  hypothetical  protein  partial  mRNA,  length  =  795 
gb  BQ790640.1  M3-l-5-_8-6-01-8-28_PM.abl  whole  midgut  #1  (WMGl),  Ae.  aegypti, 
cDNA  5',  mRNA  sequence,  length  =  493 
gb  AY432167.1,  Ae.  aegypti,  conserved  unknown  mRNA  sequence,  length  =  735 
ref  XM_001663796.1,  Ae.  aegypti,  hypothetical  protein  partial  mRNA,  length  =  1,602 
gb  DV340302.1,  Ae.  aegypti,  infected  with  Plasmodium  gallinaceum  Ae.  aegypti  cDNA 
clone  NABU360,  mRNA  sequence,  length  =  807 
gb  EF173371.1,  Ae.  aegypti,  clone  BAG  ND22N19,  complete  sequence,  length  =  146,563 
gb  BC046262.1,  X.  laevis,  HSC  protein  70,  complete  cds,  length  =  2,202 
ref  XM_001861402.1 .  Culex  pipiens  quinquefasciatus,  HSP  70  B2  partial  mRNA,  length  = 
2,072 


gb  AY431470.1,  Ae.  aegypti,  cytosolic  large  ribosomal  subunit  L19  mRNA  sequence, 
length  =  832 

ref  XM_001661916.1,  Ae.  aegypti,  ribosomal  protein  L36  pai'tial  mRNA,  length  =  490 
gb  AY432710.1,  Ae.  aegypti,  cytosolic  large  ribosomal  subunit  L27A  mRNA  sequence, 
length  =  658 

refXM_001658972.1,  Ae.  aegypti,  40S  ribosomal  protein  S3  partial  mRNA,  length  =  927 
gb  AY552053.1,  Ae.  aegypti,  40S  ribosomal  protein  S16  mRNA,  complete  cds,  length  = 
674 

refXM_001649929.1,  Ae.  aegypti,  40S  ribosomal  protein  S24  partial  mRNA,  length  =  597 

refXM_001657905.1,  Ae.  aegypti,  60S  ribosomal  protein  L31  paitial  mRNA,  length  =  632 
refXM_001647784.1,  Ae.  aegypti,  60S  ribosomal  protein  L12  partial  mRNA,  length  =  498 
refXM_001657661.1,  Ae.  aegypti,  60S  ribosomal  protein  L8  partial  mRNA,  length  = 

1,087 

refXM_001651406.1,  Ae.  aegypti,  60S  ribosomal  protein  L26  partial  mRNA,  length  =  704 
refXM_001659089.1,  Ae.  aegypti,  60S  ribosomal  protein  L24  partial  mRNA,  length  =  682 
refXM_001655416.1,  Ae.  aegypti,  60S  ribosomal  protein  LIO  pai*tial  mRNA,  length  =  959 

refXM_001651695.1,  Ae.  aegypti,  actin  partial  mRNA,  length  =  1,748 


nt:  738-1 
nt:  1933-2498 
nt:  1935-2501 
nt:  1935-2501, 
nt:  1551-1215 
nt:  49-139 
nt:  34-613 
nt:  1126-1415, 
nt:  613-34 

nt:  367-7 
nt:  7-367 
nt:  1,074-850 
nt:  2,391-2,225 

nt:  44-284 

nt:  571-486 
nt:  1,485-1,955 

nt:  655-861 
nt:  1,950-2,084 
nt:  2,004-2,370 

nt:  57-521 

nt:  1.149-1,198 
nt:  263-13 
nt:  5-315 

nt:  9-245 
nt:  187-444 
nt:  566-214 
nt:  214-566 
nt:  64,254-63,965 
nt:  1,258-1,502 
nt:  1,915-1,565 
nt:  1.475-647 


nt:  1-251 
nt:  300-39 
nt:  110-11 
nt:  19-299 


nt:  83-354 
nt:  471-329 

nt:  444-28 

nt:  21-385 
nt:  1-387 
nt:  956-707 

nt:  297-22 
nt:  72-242 
nt:  322-641 

nt:  18-286 
nt:  296-73 


lational  inhibitor  (Table  1).  We  also  found  19  clones 
encoding  approximately  nine  unknown  genes.  Com¬ 
parison  of  blastn  nucletide  sequences  of  clones  3  and 
32  with  Xenopus  laevis  Daudin  showed  a  similar  pu¬ 
tative  HSC  protein  70  with  76%  identity  at  the  nucle¬ 
otide  level;  however,  comparison  of  deduced  amino 


acid  sequences  of  clones  3  and  32  with  Anopheles 
gambiae  Giles  and  Culex  pipiens  Linnaeus  HSP  70 
indicated  95%  identity  at  the  protein  level.  Genes  of 
clones  4,  29,  and  66  had  80%  DNA  sequence  similarity 
to  the  protein-coding  regions  of  the  Cx.  pipiens  quin- 
quefasciatiis  Say  HSP  70  B2  partial  mRNA.  Therefore, 
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HSP  related  clones  Unknown  gene  clones 


Clones’ No.  9  17  14  16  20  36  45  49 

Non-HSP  related  clones 


Clones' No.  1  27  43  2  30  67  5  8  11  24  33  40 

Ribosomal  gene  clones 


Fig.  1.  Quantitative  RT-PCR  results  showing  the  relative  ratio  of  different  gene  expressed  (up-regulated  times)  after  42°C 
heat  shock  treatment  for  1  h  compared  with  the  23°C  control.  (A)  Different  HSP  genes  differentially  expressed  for  after  42°C 
heat  shock  treatment  1  h  compared  with  the  23°C  control.  (B)  Different  non-HSP  genes  differentially  expressed  after 
42°C  heat  shock  treatment  for  1  h  compared  with  the  23°C  control.  (C)  Different  unknown  genes  differentially  expressed 
after  42°C  heat  shock  treatment  for  1  h  compared  with  the  23°C  control.  (D)  Different  ribosomal  genes  differentially 
expressed  after  40°C  heat  shock  treatment  for  1  h  compared  with  the  23°C  control. 


clones  3  and  32,  combined  with  clones  4,  29,  and  66, 
might  be  new  HSC70  cognates  or  new  HSP70  that 
have  not  yet  been  identified  in  Ae.  aegypti.  In  addition, 
we  also  found  15  clones  that  encoded  for  Ae.  aegypti, 
cytosolic  large  ribosomal  subunit,  40S  ribosomal  pro¬ 
tein,  and  608  ribosomal  protein,  as  well  as  cytoskeleton 
protein  actin  (Table  1). 

High  Temperature  Effects  on  Relative  RNA  Expres¬ 
sion  Levels  of  Different  Genes  in  Adult  Ae.  aegypti.  To 
understand  whether  HSP  and  other  genes  in  Ae.  ae¬ 
gypti  can  be  triggered  by  exposure  to  high  environ¬ 
mental  temperature  conditions,  quantitative  RT-PGR 
analyses  were  carried  out  to  further  confirm  differ¬ 
ential  gene  expression  of  these  genes  in  response  to 
heat  shock  in  Ae.  aegypti.  Seven  HSP  and  HSG  genes 
were  identified  in  the  subtractive  library  (Table  1). 
According  to  our  quantitative  RT-PGR  data,  all  HSP 
genes  and  HSG  protein  genes  found  in  the  subtracted 
cDNA  library  were  up-regulated  (Pig.  IB) .  Gompared 
with  the  control,  the  RNA  relative  gene  expression 
level  of  HSP  family  increases  dramatically  after  42°G 
treatment  of  Ae.  aegypti  for  1  h  (Fig.  lA) .  For  example, 
expression  of  small  HSP  genes  (26-kDa  inRNA,  clone 
51 )  were  the  most  up-regulated  genes  after  1  h  of  42°G 
treatment  of  female  Ae.  aegypti,  more  than  a  39-fold 
relative  increase  over  that  found  in  the  untreated 
control  (1  h  at  23°G)  Ae.  aegypti  female  (Fig.  lA).  A 
five-  to  six-fold  increased  expression  of  HSG  70  mRNA 
(clones  22  and  44)  was  found  after  heat  shock  at  1  h 
of  42°G  treatment  of  Ae.  aegypti  compared  with  the 
untreated  control  (Fig.  lA;  Table  1). 


In  the  non- heat  shock  group,  quantitative  RT-PGR 
was  used  to  further  confirm  those  genes  were  either 
up-regulated  or  down-regulated  (Fig.  IB).  Six  known 
genes  were  validated  as  up-regulated  and  two  known 
genes  as  down-regulated  in  the  subtractive  hybridiza¬ 
tion  library  (Fig.  IB). 

For  the  unknown  genes,  the  quantitative  RT-PGR 
data  also  confirmed  some  up-regulated  or  down-reg¬ 
ulated  genes  (Fig.  IG).  Glone  47,  whose  mRNA  se¬ 
quence  was  similar  to  Ae.  aegypti  infected  with  Plas¬ 
modium  gallinaceum  Emile  Brumpt  clone  NABU360, 
showed  more  than  a  30-fold  relative  increase  when  Ae. 
aegypti  were  exposed  1  h  at  42°G  than  that  found  after 
being  held  continuously  at  23°G  (Fig.  IG).  It  is  not 
known  if  clone  NABU  360  represents  a  gene  from  Ae. 
aegypti  or  P.  gallinaceum.  Glones  3  and  32,  with  «“76% 
similarity  to  the  sequence  from  the  X.  laevis  Daudin 
HSG  protein  70,  showed  greater  than  a  15-fold  relative 
increase  at  42°G  than  that  found  in  the  control 
(Fig.  IG). 

In  response  to  heat  shock,  expression  levels  of  sev¬ 
eral  ribosomal  genes  were  also  significantly  modu¬ 
lated,  being  more  than  two-fold  up-regulated.  Quan¬ 
titative  RT-PGR  analysis  further  substantiated  the 
differential  expression  of  three  of  these  genes  in  re¬ 
sponse  to  heat  shock  treatment  (Fig.  ID).  Expression 
levels  of  Ae.  aegypti  cytosolic  large  ribosomal  subunit 
L19  mRNA  were  up-regulated  «“16-fold  compared 
with  the  non-  heat  shock  treatment  control.  Most  dif¬ 
ferent  gene  fragments  from  60S  ribosomal  protein  and 
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40S  ribosomal  protein  were  up-regulated  between 
two-  and  four-fold  compared  with  the  control. 

Discussion 

High  Temperature  Effects  on  Different  Gene  Ex¬ 
pression  in  Adult  Ae.  aegypti.  In  nature,  mosquitoes 
can  be  subjected  to  temperature  extremes  and  have 
developed  mechanisms  to  survive  these  conditions. 
We  used  suppression  subtraction  hybridization  to 
identify  potential  critical  gene  pathways  that  could  be 
targeted  as  part  of  a  control  strategy.  Quantitative 
RT-PCR  analysis  further  confirmed  the  differential 
expression  of  ==15  of  these  genes  in  response  to  heat 
shock  treatment  at  42°C  for  1  h  (Eig.  1 ) .  Twelve  known 
genes  were  validated  as  being  up-regulated  and  five 
known  genes  as  down-regulated  by  the  subtractive 
hybridization  library. 

In  general,  mosquitoes  have  a  core  group  of  genes 
that  are  activated  in  response  to  different  types  of 
stress.  Overexpression  of  HSPs  can  also  be  triggered 
by  exposure  to  different  kinds  of  environmental  stress 
conditions,  such  as  infection,  inflammation,  exposure 
of  the  cells  to  toxins,  starvation,  hypoxia,  or  water 
deprivation  (Mosser  et  al.  1988;  Yamuna  et  al.  2000; 
Roone  and  Vijayan  2002a,  b;  Spees  et  al.  2002;  Cheng 
et  al.  2003;  Sim  et  al.  200.5;  Chuang  et  al.  2007;  Sim  et 
al.  2007) .  Consequently,  HSPs  are  also  referred  to  as 
stress  proteins,  and  their  up-regulation  is  sometimes 
described  more  generally  as  part  of  the  stress  re¬ 
sponse.  In  a  previous  study,  gene  expression  of  An. 
gambiae  was  significantly  modulated  in  response  to 
O’nyong-nyong  virus  infection,  including  a  putative 
HSP  70,  HSC,  elongation  factor  la,  and  ribosomal 
protein  L35  (Sim  et  al.  2005,  Sim  et  al.  2007).  Other 
data  have  shown  an  essential  role  for  an  RNA  poly¬ 
merase  II  elongation  factor  in  the  regulation  of  heat 
shock  gene  expression  in  an  animal  model  (Gerber  et 
al.  2005).  According  to  this  study,  HSP  70,  HSC  70E, 
elongation  factor  2,  and  several  ribosomal  subunits 
(L24,  L36,  and  more)  were  also  confirmed  to  be  dif¬ 
ferentially  expressed  in  response  to  heat  shock  treat¬ 
ment. 

Increased  expression  of  small  HSPs  (sHSPs)  is 
known  to  be  a  key  regulatory  mechanism  in  extending 
tolerance  to  a  variety  of  environmental  stresses.  Eor 
example,  sHSP  are  expressed  in  ( 1 )  the  wasp  Ventiiria 
canescens  after  exposure  to  different  temperatures 
(Reineke  2005);  (2)  in  Anopheles  vectors  because  of 
interaction  with  Plasmodium  parasites  (Lefevre  et  al. 
2007),  and  (3)  Drosophila  in  response  to  injuries  and 
aging  (Morrow  et  al.  2004).  Our  quantitative  RT-PCR 
data  also  showed  significant  differences  in  the  expres¬ 
sion  of  sHSP  26  kDa  (clone  51,  «“40  times  up-regulated 
compared  with  the  control)  in  response  to  heat  shock 
treatment  of  Ae.  aegypti  female.  Eor  the  other  HSPs, 
the  relative  gene  expression  of  clones  62  and  50  (Pig. 
lA)  were  also  up-regulated  >20  times  compared  with 
the  control.  In  summary,  elevated  gene  expressions  of 
these  sHSP  and  HSPs  indicate  that  they  are  important 
genes  under  stressful  high  temperature  conditions  and 
are  potential  targets  for  RNAi  knock-down  studies. 


The  cDNAs  differentially  expressed  by  heat  shock 
conditions  were  identified  using  subtractive  hybrid¬ 
ization  library  analysis.  The  quantitative  RT-PCR  anal¬ 
yses  confirmed  that  genes  homologous  to  serine-type 
endopeptidase  (clones  9  and  26),  cytochrome  c  oxi¬ 
dase  (subunit  VIA)  (clone  16),  adp/ATP  carrier  pro¬ 
tein  (clone  20),  elongation  factor  2,  ATP  synthase 
delta  chain,  thioredoxin  reductase,  methionine  adeno- 
syl  transferase  and  translational  inhibitor,  and  some 
unknown  genes  were  abundant  in  Ae.  aegypti  female 
after  heat  shock  treatment.  Conversely,  glutamyl 
aminopeptidase  and  casein  kinase  were  down-regu¬ 
lated  in  Ae.  aegypti  female  after  heat  shock  treatment. 
Some  of  the  up- /down-regulated  genes  identified 
(HSP  gene,  serine-type  enodpeptidase)  have  previ¬ 
ously  been  shown  to  play  critical  roles  in  response  to 
other  types  of  stress  in  both  invertebrates  and  verte¬ 
brates,  Eor  example,  it  was  reported  that  mitochon¬ 
drial  dysfunction  mutations  in  a  gene  (designated 
levij)  that  codes  for  subunit  Via  of  cytochrome  c  ox¬ 
idase  (COX),  and  the  data  from  levy  mutants  showed 
a  COX-mediated  pathway  in  Dro.sophila.  Disruption  of 
this  pathway  leads  to  mitochondrial  encephalomyo- 
pathic  effects  including  neurodegeneration,  motor 
dysfunction,  and  premature  death  (Liu  et  al.  2007). 
Therefore,  relative  RNA  expressions  of  serine-type 
endopeptidase,  cytochrome  c  oxidase  (subunit  VIA), 
and  adp/  ATP  carrier  protein  were  up-regulated  more 
than  three-fold  in  response  to  heat  shock  treatment 
and  are  potential  pathways  to  target  for  knock-down 
with  RNAi. 

This  study  suggests  that  genes  expressed  in  response 
to  heat  shock  treatment  and/ or  temperature  play  an 
important  functional  role  in  Ae.  aegypti,  perhaps  en¬ 
hancing  survival  under  high  temperature  conditions. 
Identification  of  these  groups  of  gene  families  may 
provide  critical  information  needed  for  designing 
novel  control  strategies  for  medically  important  dis¬ 
ease  vectors  and  identifying  new  pathways  to  target 
for  the  development  of  genetic  molecular  pesticides 
(Pridgeon  et  al.  2008), 
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